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ABSTRACT 

Automobile radiators are heat exchangers that are 
used to transfer thermal energy from automobile engine to the 
surrounding atmosphere for the purpose of cooling the engine. 
Over 33% of heat energy generated by the engine through 
combustion is loss as heat dissipated in the atmosphere. The 
method of solution employed in this project work to solve the 
governing equations is the Galerkin-integral weighted-residual 
method, which is achieved following the steps of transforming 
the governing equations into Galerkin-integral weighted 
residual weak form, determination of interpolations functions, 
determination of element properties, assemblage of elements 
equations into domain equations and imposition of boundary 
conditions and solving of the assembled domain equations. 

The results showed that for temperature and velocity 
distributions in the radiator tubes and inlet hose to radiator as 
the number of elements is increased the more the finite 
element solution approximates the analytical solutions. 
Temperature values are observed to decrease, with increase in 
length, from 150°C to 80°C in the radiator tubes for finite 
element analysis, analytical, and ANSYS software used; and 
the finite element solutions exactly approximate analytical 
solutions at the nodes and agree with the ANSYS result. For 
velocity distribution in the radiator tube diameter, at the tube 
walls the no-slip boundary conditions are satisfied with 
velocity increasing from the wall at velocity of 0 to the 
midsection at velocity of 50.195m/s; while for the inlet hose 
diameter, velocity increases from wall at velocity 0 to the 
maximum at the midsection velocity 669.269m/s. Finally, the 
finite element analysis method can be used to determine how 
temperature will be distributed during radiator design stage in 
order to improve on its efficiency. 

Keywords— Radiator Tube, Finite Element, Temperature 
Distribution, Velocity Distribution 


I. INTRODUCTION 

Automobile engines are internal combustion 
engines that generate huge amount of heat when air-fuel 


mixture is combusted in the combustion chamber. 
Temperatures of metal components such valves, cylinder, 
piston and cylinder head around the combustion chamber 
can exceed 538°C(Pathade et al., 2017). Therefore, 
automobile radiators, which are heat exchangers that are 
used to transfer thermal energy from automobile engine to 
the surrounding atmosphere for the purpose of cooling the 
engine. Over 33% of heat energy generated by the engine 
through combustion is loss as heat dissipated in the 
atmosphere (Sathyan, 2016).According to Gangireddy and 
Kishore (2017), to increase the surface area available for a 
radiator with its surroundings, multiple fins are usually 
attached in contact with the radiator tubes through which 
pumped liquid flow through. Air or other exterior fluid in 
contact with the fins carries off heat. More so that radiators 
used for vehicle engine cooling are either down-flow in 
which the direction of coolant flow is vertical or cross-flow 
with horizontal flow. In their review, Patel and Dinesan 
(2014)in their parametric study used CFD to compare heat 
transfer and pressure drop of radiators with different 
parameters for optimum performance. The results showed 
that CFD results have high correlation level with the actual 
experimental results. In the work of Priyadharshini (2016), 
Pro-E 3-D modeling and finite element analysis were used 
in analyzing a radiator. The result showed Alalloy 
temperature of variation above 145°C and heat flux value of 
0.142W/m 2 ;for increasing length of fins, heat transfer value 
will be in the range of 740 watts while for copper alloy 
temperature was in the range of 145°C with heat transfer 
range above 1000 watts. And 150°C for brass with heat 
transfer in the range above 635watts. 

Ng et al. (2005) noted that for the purposes of 
increasing the boiling point of the coolant and preventing 
corrosion in the cooling systems, coolants are generally a 
mixture of water, ethylene glycol (anti-freezing agent), and 
possibly various corrosion inhibitors. It is noted that the use 
of glycol mixture generally reduces the heat transfer 
performance compared with pure water. Oliet et al. (2007) 
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observed that the heat transfer and the performance of a 
radiator are parameters that strongly affected by air and 
coolant mass flow rate; that as air and coolant flow 
increases cooling capacity also increases; and that when air 
inlet temperature increases, heat transfer and cooling 
capacity decreases. According to Romanov and Khozeniuk 
(2016), thermal load on an engine can be reduced by 
controlling the convective component of heat flows using 
the changing flow pattern of the coolant in the cooling 
system. 

The dynamics of flow in a radiator are known to 
be governed by sets of partial differential equations (PDEs) 
along with boundary conditions which are usually 
formulated to simplify the PDEs. By similarity 
transformations, the PDEs are transformed into ordinary 
differential equations (ODEs) which can then be solved 
numerically when the boundary conditions are applied 
(Reddy, 1998). 

One major problem of radiator is low rate of heat 
transfer which usually results in low efficiency challenges. 
This low heat flow rate depends on the distribution of 
temperature and the manner of velocity distributions in the 
radiator. These have prompted a lot of numerical models 
being developed to predict the cooling air flow, either in 2D 
or 3D using off-the-shelf commercial Computational Fluid 
Dynamics (CFD) software like Fluent, Vectis, StarCD and 


StarCCM+. In all the literature reviewed mathematical 
model for temperature and velocity distributions were not 
developed using finite element analysis (FEA) 
method(Pang et al., 2012; Pathade et al., 2017; Gangireddy 
and Kishore, 2017: Ng et al.,2005). Therefore in this study, 
finite element analysis method is used to develop 
mathematical models that will predict temperature and 
velocity distributions in a car radiator. 

II. MODEL GEOMETRY 

In this study 1 -dimensional flow at inlet and outlet 
hoses of the radiator as well as through the radiator tubes 
from top to bottom, down-flow radiator are considered (Fig. 
1 ). 

The fluid is assumed to enter the tube at inlet (x, 
yd with uniform inlet velocity (ui = 0), where it comes in 
contact with the inlet tube surface with viscous reaction 
between fluid and the tube thus assumed to set up a one 
dimensional (1-D), steady, uniform inward flow. The 
viscous force between the tube surface and the adjacent 
fluid layer tends to slow down the fluid velocity, which 
results in developing velocity gradient in the flow field. As 
the fluid is set in motion, pressure differential (drop) is also 
set up in the direction of flow. 



Fig. 1: Radiator tube geometry 
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Fig. 2 shows the flow characteristic of an ideal 
fully-developed viscous flow profile in the flexible hose 
tubes and radiator tubes. As the fluid enters at inlet XiOr y{ 


as the case may be, it sustains a uniform flow at this point 
but becomes laminar and fully-developed significant 
distance from the inlet with parabolic flow profile. 



Fig. 2: radiator tube and flexible tube domain flow 


2.1 Domain Discretization 

In this study, uniform and linear elements were 
chosen for the domain discretization. The radiator tube 
domain Q!(-l<x<l)is subdivided into N number of 
linear elements mesh along the y axis; while the flexible 


Radial heat conduction 



tube domain fl 2 (-1 < x < 1) is also subdivided into N 
number of linear elements mesh along the y axis. The 
radiator tube flow field and wall are subdivided into N = 4, 
linear and uniform elements mesh as in Fig. 3(a) and (b). 



Thickness 


(a) 


(b) 


Fig. 3: (a) Radiator tube flow field (b) Radiator thickness meshes 

III. MATHEMATICAL For this study, flow of fluid in the flexible tube 

FORMULATIONS and ra( Fator tube cores are governed by the continuity, 

Navier-Stokes and energy partial differential equations. 

3.1 Governing Equations 
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Continuity equation: 

du du du 

— + — + — = 0 ( 1 ) 
dx dy dz 


x-momentum equation: 
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dt dx dy dz 


dp 

■—-+ U 

dx 


f ^2 r{2 rsl \ 

8 u 8 u 8 u 




■ + ■ 


dy 1 dz 1 


+ PS X 
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Energy equation 


pc v 


dT dT dT 

U -bV- + W - 

dx dy dz 


= k 


d 


dx 


dT_ 

dx 


+ 


err crr_ 

dy 2 dz 1 


+ <!> 


(5) 


3.2 Relevant Assumptions 

The following assumptions are applied to simplify 
continuity, Navier-Stokes and equations: 

(i) One dimensional in the hose and radiator (the 
flow in the radiator is in y-direction with the 
flow in x-axis direction assumed 
insignificantly negligible), 

(ii) Incompressible (i.e. density is assumed 
constant), 

(iii) Steady, 

(iv) Laminar and fully-developed, 

(v) Newtonian and isotropic fluid flow, 

(vi) Body forces (gravitational and inertia ) are 

negligible, 

(vii) No-slip condition exists at tube inner surfaces, 
(viii) Elements are linear and uniformly spaced 

between nodes, 


(ix) Heat conduction is steady with no heat 
generation, 

(x) Heat conduction is one dimensional, along the 
y-axis only, 

(xi) Thermal conductivity is constant (isotropic 
material), 

(xii) Element length is uniform over domain region, 

(xiii) Heat transfer by radiation negligible 

By modeling the radiator this way is to enable us to 
know how temperature varies in the direction of heat 
transfer and to develop a simple approximate mathematical 
model to predict the temperature distribution in one 
dimension in the radiator. 

3.3 Initial inlet and Boundary Conditions 

The following initial inlet and essential boundary 
conditions are specified at the tube-fluid and radiator-fluid 
boundary interfaces as well as at inlet and outlet as follows: 


Boundary conditions at flexible hose tube inlet 

v(y = 1) = v(_y = —1) = 0;V ave = 0.20 —= 18-^- 

v in 


Boundary conditions at flexible hose tube outlet 

V(x,y i ) = 0.018 — -,p A =23^ 

.v in 
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Boundary conditions at radiator tube inlet and outlet 


u(x = 1) = u(x = —1) = 0; 
T a =80 °C,T l = 150° C; 
T a =30 °C,T S =150 °C 


By applying the assumptions and boundary conditions to equations (1), (2), (3), (4) and (5) reduce them to the followings 
governing equations (6), (7) and (8): 


x-momentum equation applicable to tubes: 

_ 1 dp u( d 2 u" 

0 =-- + — — T 

p dx p\dy J 


( 6 ) 


y-momentum equation applicable to radiator: 


dp 

0 —-h U 

dy 


r ii \ 

d v 


dx 2 


Energy equation 
d 2 T 


0 = Aq + k 


dy 2 


(7) 


( 8 ) 


IV. METHOD OF SOLUTION 

The method of solution employed to solve the 
close-formed governing equations (6), (7) and (8) is 
Galerkin-integral weighted-residual method. 


d f, dr) 

— kA — 
dy\ dy J 


+ qA = 0 


4.1 Heat transfer in Wall of Length L 

For steady, 1-dimensional heat conduction along 
wall, length L (y-direction), the above equation reduces to: 


(9) 


Subject to the following boundary conditions: 

m /y 

T(y = 0) = 150° C\T{y = L) = 80°C;VA = 0.20 — ; p y = 18— (10) 

s m 


Finite element models for 3,4 and 6 elements mesh are (see Appendix A): 


For 3 elements mesh: 

T x =150 °C\T 2 =12 TC,T 3 =103 °C,T 4 =80 °C 
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For 4 elements mesh: 

T x =150 °C\T 2 - 134°C;T 3 =117 ° C\T 4 =99 °C\T 5 =80°C 

For 5 elements mesh: 

T x =150 °C;T 2 =138 ° C;T 3 = 126° C;T 4 =112 ° C\T 5 = 96° C;T 6 =80 °C 

4.2 Heat Transfer in Wall of Thickness S 

For steady, one dimensional heat conduction across the wall thickness S, x-direction, the governing equation (3.13) 
becomes: 


dx 


kA— 

dx 


+ qA — 0 


Subject to the following boundary conditions: 


m A J 

T(x = 0) = 150° C;T(x = S) = 30° C;V = 0.018-; p x = 23 — 

s m 

Finite element models for 3,4 and 5 elements mesh are (see Appendix A): 

For 3 elements mesh: 


T x =150"C;r 2 =110° C;r 3 =10°OJ 4 =30 °C 

For 4 elements mesh: 


T x = 150° C\T 2 =121 °CJ 3 = 92° C\T 4 =61° C\T 5 =30 °C 

For 5 elements mesh: 

T x =150 °C;T 2 =128 ° C;T 3 =106° C;T 4 = 82° C;T 5 = 56° C;T 6 =30 °C 


( 12 ) 


4.3 Governing Equation for Flow in Radiator Tube 

The governing equation for inlet tube flow in one-dimension is: 


dx 2 jU 

subject to the following initial inlet and boundary conditions: 
v(y = -1) = 0;v(y = +1) = 0;V X = 0.20 —; p = 18 A 

V 111 

Finite element models for 3and 4 elements mesh are (see Appendix A): 


(13) 
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For 3 elements: 



f Ph 2 ) 



f PJr ) 


II 

o 

II 

to 

2 -“ j 

;v 3 

= 2 

J 

o 

II 
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4.4 Governing Equation for Flexible Hose at Radiator Inlet 

The governing equation for radiator inlet tube in one-dimension is: 

d 2 u P x 
— T -— = 0 
dy ju 

subject to the following initial inlet and boundary conditions: 

u(y =+l)=0-My = -V>=0-,V x = 0.018 — ; p = 23 ^- 

s m 

Finite element models for 3and 4 elements mesh are (see Appendix A): 


(14) 


For 3 elements 


u x — 0; u 2 =2 

For 4 elements 
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V. RESULTS AND DISCUSSION 

The results and discussion presented here are for 
temperature distributions along radiator tube length (L) and 
thickness (S) as well as velocity distributions along radiator 
tube length (L) and in the radiator inlet hose. 

For temperature distribution in the radiator tube 
length, Fig. 4 is a graph of temperature against tube length 
comparing finite element solutions for 3 (diamond), 4 
(square) and 5 (triangle) elements with analytical solution 
(cross); the graph reveals that the more the number of 
elements is increased the more the finite element solution 
approximates the analytical solutions. Fig. 5 shows 
temperature distribution long tube length using ANSYS 
R16.2 software. Temperature is seen to decrease from 


150°C (red portion) to approximately 80°C (blue portion), 
which agrees with the result obtained by analysis. Fig. 6is a 
graph of temperature against tube length comparing finite 
element temperature distribution pattern (diamond) with 
heat flow distribution pattern (triangle) along the tube 
length respectively. In Fig. 6, temperature is seen to 
decrease with increase in length while heat flow increases 
with increase in length. For temperature distribution along 
the tube thickness, Fig. 7 is a graph of temperature against 
tube thickness comparing finite element solutions for 3 
(diamond), 4 (square) and 5 (triangle) elements with 
analytical solution (cross); the graph reveals that the more 
the number of elements is increased the more the finite 
element solution approximates the exact solutions. Fig. 8is 
a graph of temperature against tube thickness comparing 
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finite element temperature distribution pattern (diamond) 
with heat flow distribution pattern (triangle) along the tube 
thickness respectively with temperature shown to decrease 
with increase in length while heat flow increases with 
increase in length. 

For velocity distribution in the radiator tube 
diameter, Fig. 9 is a graph of velocity against tube length 
comparing finite element solutions for 3 (square) and 4 
(diamond) elements with analytical solution (cross); the 
graph reveals that more the number of elements is increased 
the more the finite element solution approximates the exact 
solutions. From Fig. 9, at the tube walls the no-slip 


boundary conditions are satisfied with velocity increasing 
from the wall velocity = 0 to the midsection with velocity = 
50.195m/s. While for velocity distribution in inlet hose 
diameter, in Fig. 10 a graph of velocity against inlet hose 
diameter, 2 elements solutions (triangle), 3 elements 
solutions (diamond) and 4 elements solutions (circle) are 
compared with the analytical solution (square).The graph 
reveals the finite element solution approximates the 
analytical solution as the number of elements is increased 
with the no-slip boundary conditions being satisfied at the 
walls as velocity increases from 0 to the maximum at the 
midsection with velocity = 669.269m/s. 


Radiator Tube Length Finite element and 
Analytical Temperature Distribution 
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Fig. 4: Finite element and analytical temperature distribution along tube length 
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Fig. 5: Finite element ANSYS temperature distribution along tube length 
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Fig.7: Finite element and analytical temperature distribution along tube thickness 
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Fig. 10: Velocity distribution radiator inlet hose 


For radiator tube, Table 1 shows temperature 
distribution at the element nodes for 3 (3E), 4(4E) and 
5(5E) elements and the corresponding analytical (AS) 
values. In Table 1, the 5 elements nodal values approximate 
with those of the analytical nodal point values. Table 2 
shows temperature distribution for 5 (5E) elements, the 
corresponding analytical (AS) values and the heat flux 
values. Temperature distribution values are seen to decrease 
with tube length as the heat flux increases. For tube 
thickness, Table 3 indicates temperature distribution values 
at element nodes for 3 (3E), 4(4E) and 5(5E) elements and 
the corresponding analytical (AS) values. In Table 3, the 5 
elements nodal values approximate with those of the 
analytical nodal point values. Table 4 shows temperature 


distribution for 5 (5E) elements, the corresponding 

analytical (AS) values and the heat flux values. 
Temperature distribution values are seen to decrease with 
tube length as the heat flux increases. 

For radiator tube, Table 5 shows velocity 
distribution at the element nodes for 3 (3E) and 4(4E) 
elements and the corresponding analytical (AS) values. In 
Table 5, the 4 elements nodal values approximate with 
those of the analytical nodal point values. For inlet hose, 
Table 6 shows velocity distribution at the element nodes for 
3 (3E) and 4(4E) elements and the corresponding analytical 
(AS) values. In Table 6, the 4 elements nodal values 
approximate with those of the analytical nodal point values. 


Table 1: Radiator tube length finite element and analytical temperature nodal values 


Nodes 

3E (°C) 

4E (°C) 

5E (°C) 

AS (°C) 

0.00 

150.000 

150.000 

150. 000 

150.000 

0.20 

127.000 

134.000 

138.000 

138.488 

0.40 

103.000 

117.000 

126.000 

125.733 

0.60 

95.000 

99.000 

112.000 

111.733 

0.80 

88.000 

90.000 

96.000 

96.488 

1.00 

80.000 

80.000 

80.000 

80.000 


3E: 3 elements solution; 4E: 4 elements solution; 5E: 5 elements solution; AS: analytical solution 


Table 2: Finite element and analytical temperature distribution and heat flow values 


Nodes 

5E (°C) 

AS (°C) 

Q(J/m 2 ) 

0.00 

150. 000 

150.000 

30.545 

0.20 

138.000 

138.488 

34.035 
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0.40 

126.000 

125.733 

37.525 

0.60 

112.000 

111.733 

41.015 

0.80 

96.000 

96.488 

44.505 

1.00 

80.000 

80.000 

47.995 


5E: 5 elements solution;AS: analytical solution; Q: heat fluxes 


Table 3: Finite element and analytical temperature nodal values 


Nodes 

3E (°C) 

4E (°C) 

5E (°C) 

AS (°C) 

0.00 

150.000 

150.000 

150. 000 

150.000 

0.06 

110.000 

121.000 

128.000 

128.48 

0.12 

70.000 

92.000 

106.000 

105.733 

0.18 

57.000 

61.000 

82.000 

81.733 

0.24 

43.000 

46.000 

56.000 

56.488 

0.30 

30.000 

30.000 

30.000 

30.000 


3E: 3 elements solution; 4E: 4 elements solution; 5E: 5 elements solution; AS: analytical solution 


Table 4: Finite element, analytical temperature and heat flow nodal values 


Nodes 

5E (°C) 

AS (°C) 

Q(J) 

0.00 

150. 000 

150.000 

58.595 

0.20 

128.000 

128.488 

62.085 

0.40 

106.000 

105.733 

65.575 

0.60 

82.000 

81.733 

69.056 

0.80 

56.000 

56.488 

72.555 

1.00 

30.000 

30.000 

76.045 


5E: 5 elements solution;AS: analytical solution; Q: heat fluxes 


Table 5: Radiator tube FE and exact solutions velocity nodal values 


b 

3E 

4E 

AS 

1.000 

0 

0 

0 

0.875 

0.056 

0.113 

0.113 

0.750 

0.105 

0.210 

0.210 

0.625 

0.146 

0.293 

0.293 

0.500 

0.180 

0.361 

0.361 

0.375 

0.207 

0.413 

0.413 

0.250 

0.225 

0.451 

0.451 

0.125 

0.237 

0.473 

0.473 

0.000 

0.240 

0.481 

0.481 


Table 6: Radiator flexible hose FE and Exact velocity nodal values 


b 

3E 

4E 

AS 

1.000 

0 

0 

0 

0.875 

0.056 

0.113 

0.113 

0.750 

0.105 

0.210 

0.210 

0.625 

0.146 

0.293 

0.293 

0.500 

0.180 

0.361 

0.361 

0.375 

0.207 

0.413 

0.413 
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0.250 

0.225 

0.451 

0.451 

0.125 

0.237 

0.473 

0.473 

0.000 

0.240 

0.481 

0.481 


3E: 3 elements solution; 4E: 4 elements solution; AS: analytical solution 


VI. CONCLUSION 

From this empirical finite element analysis of 
radiator temperature and velocity distribution, it is here 
concluded that: 

• Empirical analysis on a radiator was carried out to 
obtain required parameters. 

• Finite element method was used to solve 
Continuity and Navier Stokes equation. 

• FEA was used to determine velocity and 
temperature distribution in radiator. 

• Analytical solutions were derived and compared 
with FE solutions. 

Furthermore, this study on finite element analysis of 
radiator temperature and velocity distributions revealed that 
as the number of elements in the analysis are increased FE 
results tend to converge fast to the analytical results. Also, 
the finite element nodes values are observed to be the same 
with analytical values at the nodes, which is a good 
indication that the FEA method is suitable for determination 
of temperature distribution. This is further validated by the 
ANSYS R16.2 results. Consequently, the FEA method can 
be confidently used to determine how temperature will be 
distribute during radiator design stage. 

RECOMMENDATIONS 

There is need for possible further research by 
using other methods of investigation, geometry of radiator 
and by considering the convection heat transfer by air 
current. 
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Appendix A 


(1) Transformation of governing equation into Weak form 

The derivation of element equations for typical element e in the mesh follows these steps (Reddy, 1995): 


d_ 

dy 


f 


, , dT 
kA — 

V dy 


+ qA — 0 


subject to the following boundary conditions: 


(Al) 


r(y = 0) = 150 °C;T(y = H) = 80°C 


(A2) 


The residual, R, of the governing equation (Al) is: 




, , dT 
kA — 

V dy 


+ qA — 0 


d_ 
dy 

The Galerkin-weighted residual integral equation (A4) is: 


o=l' 


[ w 

—kA 

r d 2 T^ 

i 2 

+ Aq(y) 

J 

Q 


ydy J 



dn 


The Weak form of the Galerkin-weighted residual integralequation (A5) is: 

dT 


0 = 1 


7 . dw dT 

kA -b wAq 

dy dy 


dCl — wkA - 


dy 


The finite element equation is: 

0 = 1 


= J 7 Z 1 M 


dy ji ' dy 
h dy/. dyz e . 


\ 

e dT T 

dy - 

kAy/: — 

) 

L d y \ L 


j = 1 o 


dy dy 


dy + Aq^y/Uly 


i a e dT 
kAy/i — 
dy 


0 h 


J L 


In matrix form the finite element model is written as: 
[K e ]{T e } = {f e }-{Q e } (A7) 

(2) The analytical solution are: 


T (y) = -S!L + 

y ’ 2k 



/ \ 

2 

2 _ 

y 


L 

UJ 



y 


+ (t l -t 0 )^ + t o 


Q(y) = kA 


j rAL dT(y) _ qLA 

1-2 

f y\ 

dy 2 


U )_ 


— ~(t l -t 0 )—+t o 

L y 0 L 


(A3) 


(A4) 


(A5) 


(A6) 


(A) 

(A9) 
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Appendix B 


Problem data 
Empirical values 


Observed parameters 

Symbols 

Quantity 

Unit 

Tube material 

- 

Aluminium 


Inlet temperature 

T l 

150 

°c 

Outlet temperature 

T 0 

80,30 

°c 

Specific heat 

c p 

4.187 

kJ/kgK 

Thermal conductivity 

k 

0.66 

W/mK 

Density 

P 

1000 

kg/m 3 

Dynamic viscosity 

P 

1.793 x 10' 3 

kg/ms 


Finite element analysis parameters 


Observed parameters 

Symbols 

Quantity 

Unit 

Specified pressure at tube 

Py 

0.018 

N/m 2 

Specified pressure at hose 

Px 

0.23 

N/m 2 

Inlet diameter of tube 

D, 

10 

Mm 

Outlet diameter of tube 

d 2 

11.25 

Mm 

Tube thickness 

s 

0.0625 

Mm 

Element length, tube 

h y 

0.5 

Mm 

Element length, hose 

h x 

0.005 

Mm 
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